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Resumen 


La sequía es uno de los desastres naturales más frecuentes en el centro-norte 
de México, y su identificación y caracterización son cruciales para tratar de 
mitigar sus efectos. Se analizó la evolución de las sequías meteorológicas en 
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el estado de Zacatecas, México, para el periodo de 1961 a 2012 utilizando tres 
índices multiescalares: el Índice de Precipitación Estandarizado (SPI), basado 
en datos de precipitación; el Índice de Reconocimiento de Sequía (RDI), y el 
Índice de Precipitación-Evapotranspiración Estandarizado (SPEI), basados en 
datos de precipitación y evapotranspiración. El patrón temporal de las sequías 
se analizó en tres escalas de tiempo (3, 6 y 12 meses). Los índices se 
compararon empleando los coeficientes de correlación de Pearson y de 
Spearman. Los principales periodos secos fueron 1982-1983, 1998-2000 y 
2010-2012. En la escala de tres meses, el SPI detectó el menor número de 
eventos de sequía (41), con una duración media de 2.5 meses; y en las escalas 
de 6 y 12 meses, el RDI detectó el menor número de eventos de sequía, 32 y 
16, respectivamente, con duración media de 3.0 y 6.1 meses, 
respectivamente. En la escala de tres meses, el SPEI detectó el mayor número 
de eventos extremos, y en las escalas de 6 y 12 meses fueron el RDI y el SPI, 
respectivamente. En las tres escalas de análisis, el SPEI presentó el mayor 
porcentaje de frecuencia promedio de sequías severas. Los índices SPI vs. RDI 
presentan la correlación más alta (> 0.92), y el SPI vs. SPEI la correlación 
más baja. Las áreas más afectadas por las sequías son la región noroeste, 
centro y sur de Zacatecas, donde se ubican las principales zonas agrícolas de 
temporal y de riego del estado. Los resultados encontrados pueden ayudar 
principalmente en la toma de decisiones sobre los patrones de cultivo y la 
asignación óptima de agua en el Distrito de Riego 034 Estado de Zacatecas. 


Palabras clave: sequía meteorológica, SPI, RDI, SPEI, regiones semiáridas, 
severidad y duración de sequías, evapotranspiración potencial, coeficientes de 
correlación de Pearson y de Spearman. 


Abstract 


Drought is one of the most frequent natural disasters to occur in north-central 
Mexico and its identification and characterization are crucial for mitigating its 
effects. This study analyzed the evolution of meteorological drought over the 
period 1961-2012 in the state of Zacatecas, Mexico, using three multiscale 
indices: the Standardized Precipitation Index (SPI) based on precipitation 
data, the Reconnaissance Drought Index (RDI), and the Standardized 
Precipitation-Evapotranspiration Index (SPEI) based on precipitation and 
evapotranspiration data. The temporal pattern of the drought was analyzed at 
three time scales (3, 6 and 12 months). The indices were compared using the 
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Pearson and Spearman correlation coefficients. The main dry periods were 
1982-1983, 1998-2000, and 2010-2012. At a 3-month scale, the SPI detected 
the lowest number of drought events (41), with an average duration of 2.5 
months; and at 6- and 12-month scales, the RDI detected the lowest number 
of drought events, 32 and 16, respectively, with an average duration of 3.0 
and 6.1 months, respectively. The SPEI detected the greatest number of 
extreme events at the 3-month scale, and the RDI and the SPI did so at the 
6- and 12-month scales, respectively. For the three scales of analysis, the SPEI 
had the greatest average frequency percentage of severe droughts. The SPI 
vs RDI had the highest correlation (> 0.92), and the SPI vs SPEI had lowest 
correlation. The areas that were most affected by droughts were the 
northwestern, central, and southern regions of Zacatecas, where the principal 
rainfed and irrigated agriculture areas of the state are located. The findings 
can primarily help decision-making regarding cropping patterns and the 
optimal allocation of water in Irrigation District 034 in the state of Zacatecas. 


Keywords: Meteorological drought, SPI, RDI, SPEI, semiarid regions, severity 
and duration of droughts, potential evapotranspiration, Pearson and Spearman 
correlation coefficients. 
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Introducción 


El centro y norte de México se encuentran dentro de la franja de alta presión 
de latitud norte, por lo que estas zonas son áridas y semiáridas, que coinciden 
en latitud con los grandes desiertos de África y Asia. En consecuencia, tales 
regiones no disponen de agua en cantidades suficientes y de forma histórica 
han sido afectadas por las sequías. En el siglo XX, México presentó cuatro 
grandes periodos de sequía: 1948-1954, 1960-1964, 1970-1978 y 1993-1996, 
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así como una sequía severa en 1998, las cuales afectaron en particular a los 
estados del norte del país (Cenapred, 2007). De acuerdo con Escalante y Reyes 
(1998), el estado de Zacatecas, ubicado en el centro-norte de México, está 
entre los estados con más daños a la agricultura y la ganadería por las sequías 
ocurridas durante el periodo de 1988 a 1995. 


Recientemente, en México se han registrado varias sequías de magnitud 
considerable entre los años 2000-2003, 2009 y 2011-2012. Entre los años 
2000 y 2003, 18 de los 32 estados de México presentaron sequía; los estados 
del norte, incluyendo Zacatecas, fueron los más afectados. Las pérdidas 
económicas ascendieron a más de 1 800 millones de pesos (188 millones de 
dólares), cerca de un millón de hectáreas de cultivo fueron afectadas y se 
perdieron más de 13 mil cabezas de ganado (Cenapred, 2002; Cenapred, 
2003; Cenapred, 2004). En 2009, México sufrió la segunda peor sequía en 60 
años: los daños ascendieron a 3 081 millones de pesos (229 millones de 
dólares) y la superficie de cultivos y pastizales afectada fue de 384 540 
hectáreas. Zacatecas estuvo entre los cinco estados con mayor superficie 
afectada (Semarnat, 2012). A finales de 2011 y hasta mediados de 2012, más 
de 60% de México sufrió la peor sequía en los últimos 70 años, clasificada 
según el Monitor de Sequía de América del Norte (NADM, por sus siglas en 
inglés) de intensidad severa a excepcional (Semarnat, 2012). Además, según 
el Departamento de Agricultura de los Estados Unidos de América (USDA, por 
sus siglas en inglés), esta sequía generó pérdidas económicas en el sector 
agrícola de más de 16 000 millones de pesos (1 300 millones de dólares) 
(Herron, 2013). El estado de Zacatecas fue de los más afectados por este 
desastre natural. Bajo tal contexto, es importante conocer el comportamiento 
puntual y regional de la sequía en el estado de Zacatecas para poder establecer 
medidas de mitigación. 


Las sequías pueden tener repercusiones graves y duraderas en los sistemas 
naturales y humanos, que incluyen desastres humanitarios, pérdidas 
económicas y tensiones en los ecosistemas naturales (Touma, Ashfaq, Nayak, 
Kao, 8: Diffenbaugh, 2015). La sequía es impulsada por variaciones climáticas 
y por condiciones hidrológicas adversas (Li, Liang, Yu, € Acharya, 2014). 
Aunque la sequía es la consecuencia de una reducción natural en la cantidad 
de precipitación en un periodo de tiempo, su severidad se evalúa en términos 
de su duración, época de ocurrencia, extensión territorial afectada y, 
fundamentalmente, en los impactos que causa a las actividades humanas y 
agrícolas, así como al medio ambiente (Caparrini 8 Manzella, 2009). La sequía 
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ha sido clasificada en cuatro categorías: meteorológica, hidrológica, agrícola y 
socioeconómica (Samaniego, Kumar, € Zink, 2013). 


La complejidad para definir una sequía se debe a la dificultad de cuantificar su 
magnitud o severidad, pues por lo habitual ésta se identifica por sus efectos 
sobre diferentes sistemas (recursos hídricos, agricultura, ecología, pérdidas 
económicas, etc.); sin embargo, no existe una variable física que permita 
medir en específico la severidad de la sequía. La cuantificación de la sequía 
por lo general se realiza empleando índices de sequía (Asadi-Zarch, 
Sivakumar, 8 Sharma, 2015), los cuales se basan en información de 
precipitación, temperatura, evapotranspiración, caudales en ríos u otras 
variables hidrometeorológicas medibles, y que cuantifican el riesgo de sequía. 
En la actualidad existe una gran cantidad de índices para cuantificar las 
sequías; los más utilizados a nivel mundial son el Índice de Precipitación 
Estandarizado (SPI, por sus siglas en inglés) desarrollado por McKee, Doesken 
y Kleist (1993), que consiste en un enfoque probabilístico de la precipitación, 
y el Índice de Severidad de la Sequía de Palmer (PDSI, por sus siglas en 
inglés), desarrollado por Palmer (1965), que analiza el balance de agua en el 
suelo. Recientemente, la Organización Meteorológica Mundial (WMO, por sus 
siglas en inglés) presentó el SPI como el índice de sequía universal (Hayes, 
Svoboda, Wall, £ Widhalm, 2011; WMO, 2012). La ventaja principal del SPI 
es que se puede calcular en diferentes escalas de tiempo (por lo habitual en 
1,3,6,9, 12, 24 y 48 meses) y representar las variaciones hidrológicas de los 
déficits de precipitación. Sin embargo, el SPI tiene una desventaja: sólo utiliza 
información de precipitación, sin tener en cuenta otras variables climatológicas 
que juegan un papel importante durante el desarrollo de un evento de sequía 
(Teuling et al., 2013). Recientemente, se han desarrollado índices que 
incorporan la evapotranspiración para evaluar la sequía. Por ejemplo, Tsakiris 
y Vangelis (2005), propusieron el Índice de Reconocimiento de Sequía (RDI, 
por sus siglas en inglés), que analiza el cociente entre la precipitación 
acumulada y la evapotranspiración potencial acumulada (ETP) desde un punto 
de vista probabilístico; por su parte, Vicente-Serrano, Beguería y López- 
Moreno (2010), desarrollaron el Índice de Precipitación-Evapotranspiración 
Estandarizado (SPEI, por sus siglas en inglés), que consiste en un balance 
hídrico climático mensual. 


La complejidad del fenómeno de sequía hace que ningún índice de sequía se 
pueda aplicar en todas las regiones o sistemas naturales. Por lo tanto, para un 
análisis de sequía es conveniente considerar más de un índice, a fin de 
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examinar la sensibilidad y precisión de cada uno de ellos, y explorar su 
comportamiento ante condiciones y finalidades específicas. La caracterización 
de un evento de sequía implica la identificación de la fecha de inicio y de 
finalización, su duración, magnitud (severidad) e intensidad (McKee et al., 
1993;  Kavalieratou, Karpouzos, 8  Babajimopoulos, 2012). Otras 
características de las sequías son la cobertura espacial (extensión superficial) 
y la frecuencia con que éstas ocurren (Kirono, Kent, Hennessy, € Mpelasoka, 
2011). 


Los objetivos de este trabajo fueron: (a) investigar las principales 
características de las sequías meteorológicas en el estado de Zacatecas 
durante el periodo de 1961 a 2012, utilizando los índices multiescalares de 
sequía SPI, RDI y SPEI; y (b) comparar el comportamiento temporal y espacial 
de los índices de sequía utilizados. 


Materiales y métodos 


Características del área de estudio 


El estado de Zacatecas está ubicado entre las coordenadas 21% 01' 00” y 25% 
07' 00” de latitud norte, y 100% 43" 00” y 104* 22' 00” de longitud oeste, y 
cuenta con una superficie de 75 284 km?. Los principales tipos de clima que 
presenta el estado son seco y semiseco (73%), y templado subhúmedo (17%). 
La temperatura media anual es de 17 *C, y las temperaturas máxima y mínima 
promedio anuales son de 30 *C (mayo) y 3 *C (enero), respectivamente. La 
precipitación media anual en Zacatecas es de 510 mm, con precipitaciones de 
300 mm en el norte y de 860 mm en el sur; las lluvias se concentran en verano 
(75%), de junio a septiembre. El clima del estado de Zacatecas es una 
limitante para la agricultura, principal actividad económica, que destina 
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alrededor de 1.74 millones de hectáreas, de las cuales 88% es de temporal 
(secano). 


Información climatológica 


Para identificar los eventos de sequía, así como para evaluar sus duraciones, 
magnitudes, e intensidades y frecuencias, en esta investigación se utilizaron 
las series de datos de precipitación y temperatura diaria de nueve estaciones 
climatológicas en el estado de Zacatecas (Figura 1), tomadas de la base de 
datos del CLICOM (Conagua, 2014), para el periodo del 1% de enero de 1961 
al 31 de diciembre de 2012. Los descriptores geográficos (coordenadas y 
altitud) de las estaciones climatológicas y los valores medios anuales 
correspondientes a precipitación, P, y temperatura, 7, se presentan en la Tabla 
1. 
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Figura 1. Estaciones climatológicas seleccionadas para el análisis de la 
sequía. 


Tabla 1. Estaciones climatológicas y su precipitación y temperatura media 
anual. 
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Latitud [Longitud | ajitud iaa | Tempe- 
Clave | Nombre msnm) Periodo | *9cI10n ratura 
(N) (0) (mm) (9C) 

: 1961- 

32052 |San Rafael 240 35” |1020 06 2 014 2012 262.8 14.4 
1961- 

32028 |Juan Aldama | 24% 16” |1030 23 1 999 2012 433.0 18.0 
1961- 

32001 | Agua Nueva |23% 47" |1020 09 1 946 2012 349.1 17.7 
Sombrerete 1961- 

2054 230 il Y d y ] 

3205 (DGE) 3038” |1030 38” |2 300 2012 551.4 16.4 
1961- 

32018 |El Sauz 230 16 |1030 06” |2 096 2012 414.7 16.0 
1961- 

32003 |Calera 22054” [1020 39 |2 097 2012 426.5 15.6 
! A y 1961- 

32030 |La Florida 22041” |1030 36 1870 2012 590.3 16.4 
Guadalupe e EAS 1971- 

32024 Victoria 22023" |1010 49 2 132 2012 374.1 16.4 
An A , 1961- 

32032 |La Villita 21036" |1030 20 1 786 2012 782.3 20.4 

Metodologia 


En esta investigación se emplearon tres de los índices de sequía más usados 
y conocidos a nivel internacional: el Índice de Precipitación Estandarizado, el 
Índice de Reconocimiento de Sequías y el Índice de Precipitación- 
Evapotranspiración Estandarizado. La caracterización de las sequías se llevó a 
cabo en seis etapas: 1) selección y control de calidad de la información 
climatológica; 2) selección de la escala de análisis; 3) cálculo de la 
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evapotranspiración potencial (ETP); 4) cálculo de los índices SPI, RDI y SPEI 
en las escalas de tiempo de 3, 6 y 12 meses; 5) análisis de las características 
de las sequías; y 6) comparación de los índices de sequía meteorológica. En la 
Figura 2 se presenta de modo esquemático el desarrollo metodológico de la 
investigación. 
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Figura 2. Diagrama esquemático del desarrollo metodológico para el análisis 
de las sequías. 
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Selección y control de calidad de la información 
climatológica 


Una etapa fundamental en el análisis de series de tiempo hidroclimatológicas, 
como es el caso del análisis de sequías, es la selección y control de la calidad 
de la información. La selección de la información climatológica empleada se 
basó en particular en cuatro criterios: (1) que la estación climatológica 
estuviera activa para poder evaluar las sequías de los últimos años; (2) que la 
longitud de los registros de precipitación y temperatura de cada estación 
climatológica fuera de al menos 40 años; (3) que el porcentaje de datos 
faltantes fuera menor o igual al 15%, para no perder representatividad en los 
datos; y (4) que las estaciones estuvieran distribuidas en toda la zona de 
estudio, a fin de tener una buena representación espacial de la sequía. 


Para el análisis de la sequía se llevaron a cabo controles de calidad preliminares 
de las series de datos de precipitación y temperatura, revisando la 
homogeneidad de las series de información climatológica y completando los 
datos faltantes. Para revisar la homogeneidad de los datos de precipitación y 
temperatura se aplicaron pruebas estadísticas paramétricas y no paramétricas 
a las series anuales. Las pruebas seleccionadas fueron la t de Student, la 
prueba de Cramer (Mirza, 1997), y las Pruebas de Homogeneidad Normal 
Estándar (SNHT, por sus siglas en inglés) y de Buishand (Hánsel, Medeiros, 
Matschullat, Petta, 8 De Mendonca-Silva, 2016). 


Para seleccionar el método de complementación de datos faltantes se probaron 
varios métodos: interpolación polinomial, promedio aritmético, inverso de la 
distancia y relación normal (Tabios III 8 Salas, 1985; Dingman, 2002); el 
método que presentó mejores resultados fue el de la relación normal (ecuación 


(1): 
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donde f,, es el dato faltante a estimar para el día, mes o año en la estación de 
análisis m; g, las estaciones vecinas; n,, el número de estaciones vecinas 
consideradas en el análisis; Pp, y Pz son la precipitación media anual en la 
estación m con información faltante y en las estaciones vecinas g, 
respectivamente, y p, es la precipitación regitrada en las estaciones vecinas el 
día, mes o año en el que falta el dato en la estación m. Este método fue 
aplicado de manera análoga a los datos de temperatura. 


Selección de las escalas de análisis de la sequía (k) 


Para investigar las características de la sequía, la escala de tiempo k juega un 
papel importante. Las escalas temporales k cortas son relevantes por ser 
representativas de sequías meteorológicas y agrícolas. Una escala temporal 
de tres meses da una estimación estacional de la precipitación, que refleja las 
condiciones de humedad a corto y mediano plazos. Una escala de seis meses 
estima la precipitación a mediano plazo, la cual puede estar asociada con 
niveles en embalses y caudales anómalos, mientras una escala temporal k más 
larga de 12 meses refleja los patrones de precipitación a largo plazo, que 
probablemente están relacionados con caudales, niveles de embalse e incluso 
niveles de agua subterránea, y permiten identificar la sequía hidrológica. 


En esta investigación fueron seleccionadas las escalas k de 3, 6 y 12 meses, 
para los índices SPI, RDI y SPEI (SPI-3, SPI-6, SPI-12, RDI-3, etc.), con el fin 
de obtener perspectivas a corto, mediano y relativamente largo plazos sobre 
las características de la sequía. 


Cálculo de la evapotranspiración potencial (ETP) 
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Para el cálculo de los índices RDI y SPEI, la evapotranspiración potencial fue 
estimada mediante el método de Thornthwaite (1948), como lo proponen 
Tsakiris y Vangelis (2005), y Vicente-Serrano et al. (2010). Este método es 
uno de los más sencillos para calcular la ETP (ecuación (2)), ya que sólo 
requiere de datos de temperatura media mensual y la latitud del lugar donde 


se desea estimar la ETP: 


mn (2) 


donde ETP es la evapotranspiración potencial mensual (mm/mes); T, 


la 


temperatura media mensual (*C), e / es el índice de calor (ecuación (3)), el 


cual se calcula como la suma de 12 valores de índices mensuales: 
FS 2) 
; (3) 


m es un coeficiente que depende del índice de calor /. 
m=6.75x 1071? — 7.71 x 10791? + 1.79 x 107?] + 0.492 (4) 


K es un coeficiente de corrección calculado como una función de la latitud 


y el mes: 


e Ea peta (5) 


donde NDM es el número de días del mes y N el número máximo de horas 


de sol (ecuación (6)). 
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n= (—)w (6) 
y w, es el ángulo horario de la salida del sol (ecuación (7)): 
w, = arccos(—tangtanó) (7) 


en la cual, y es la latitud (radianes) de la estación climatológica y 8 la 
declinación solar (radianes) calculada como: 


211] 
g = 0.4093 sen (5 = 1.405) (8) 


donde J es el día Juliano promedio del mes. 


El Índice de Precipitación Estandarizado (SPI) 


El SPI considera la precipitación como variable única de análisis y estima si en 
una región, en un determinado periodo, hay déficit o exceso de precipitación 
respecto a las condiciones normales (Hayes, Svoboda, Wilhite, £ Vanyarkho, 
1999). El cálculo del SPI se basa en registros de precipitación de largo plazo; 
Guttman (1994) señala que se necesitan de preferencia registros de 50 a 60 
años o más, aunque otros investigadores (Wu, Hayes, Wilhite, 8 Svoboda, 
2005; Li et a/., 2014) recomiendan tener al menos 30 años. El SPI se calcula 
ajustando una función de distribución de probabilidad (fdp) apropiada, 
generalmente la función de distribución Gamma, a los datos de precipitación 
mensual acumulada, para cada escala de tiempo k de interés y para cada sitio 
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de interés en el espacio; los valores del SPI se obtienen transformando los 
valores de la distribución Gamma en valores de la variable normal estándar. 
La función de densidad de probabilidad Gamma está definida por: 


g(x)= a-1e=*/boarax>0 (9) 


O 


donde a es el parámetro de forma; b, el parámetro de escala; x, la 
precipitación mensual acumulada (x > 0), en mm, y T(a) es la función Gamma. 
Los parámetros a y b pueden estimarse por varios métodos, cuyos valores por 
el método de máxima verosimilitud son (Haan, 1977): 


1 44 (10) 
EA 1+ 1+3 


Xx 
b=- CL) 
donde A es una variable auxiliar que está definida por: 
A=in(x) -- > (12) 
= In(x E Xi 
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donde n' es el número de datos de precipitación x, diferentes de cero, y x es 
el valor medio de los datos x;. 


Integrando la función de densidad Gamma (ecuación (9)) con respecto a x, se 
obtienen la función de distribución acumulada G(x): 


E Sy 
G(x) = A) e lb dx (13) 


Considerando t = x/b, entonces la ecuación (13) se reduce a la función Gamma 
incompleta: 


1 x 
G(x) = ña) pote stdt (14) 


Dado que la función Gamma no está definida para x=0, y las series de 
precipitación pueden tener valores de cero, entonces la probabilidad 
acumulada puede calcularse como  (Angelidis, Maris, Kotsovinos, € 
Hrissanthou, 2012): 


H(x) = q + (1-— q)G(x) (15) 


donde H(x) es una función de probabilidad acumulada; q, la probabilidad de 
tener valores de precipitación igual a cero, y G(x) es la probabilidad acumulada 
de la función Gamma incompleta. 


Si m es el número de ceros en la serie de precipitación, entonces q = m/n, 
donde n es el número total de datos de la serie de precipitación, según la 
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escala de análisis k. Cuando la serie no tiene valores de precipitación igual a 
cero, q =0, por lo cual, H(x) = G(x), y cuando x = 0, H(0) = q. 


Finalmente, para definir el valor del SPI, la distribución de probabilidad 
acumulada H(x) se transforma en la variable normal estándar Z de media cero 
y varianza uno (Edwards 8 McKee, 1997), utilizando la aproximación 
propuesta por (Zelen 8, Severo, 1965): 


lio : Chata t? pe 1 1 4 
o y ETRE A (160) para (16) 


<H(x) < 0.5 


O A A 
2=sPI= +0 A E ),: 


—1+d,t+dat? + dgt3 

ñ (17) 

= |In (5) para 0.5 < H(x)<1.0 
(1 HGo0) 


donde c, = 2.515517, c, = 0.802853, c, = 0.010328, d, = 1.432788, d, = 0.189269 y 
dz = 0.001308. 


Un evento de sequía ocurre cuando a cualquier escala de tiempo, el SPI es 
continuamente negativo y alcanza un valor de -1.0 o menor, y termina cuando 
el SPI llega a ser positivo (McKee et a/., 1993; Vrochidou y Tsanis, 2012). Los 
valores del SPI se pueden comparar con facilidad y de forma simultánea en los 
ámbitos espacial y temporal (Lopez-Bustins, Pascual, Pla, 8 Retana, 2013). 
Con base en el rango de valores del SPI, los eventos de sequía se clasifican de 
ligeros a extremos (Tabla 2). 


Tabla 2. Clasificación del SPI (McKee et a/., 1993). 


Valor SPI Clasificación 
> 2.00 Extremadamente húmedo 
1.50 a 1.99 Severamente húmedo 
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Valor SPI Clasificación 


1.00 a 1.49 Moderadamente húmedo 
Ligeramente húmedo (cercano a lo 

0a0.99 normal) 

O a -0.99 Sequía ligera (cercano a lo normal) 


-1.00 a -1.49 Sequía moderada 


-1.50 a -1.99 Sequía severa 


< -2.00 Sequía extrema 


El Índice de Reconocimiento de Sequías (RDI) 


El RDI se puede calcular en diferentes escalas de tiempo y su interpretación 
es similar a la del SPI (Tabla 2) (Vicente-Serrano, Van der Schrier, Beguería, 
Azorin-Molina, 8. López-Moreno, 2015; Xu et al., 2015). Entre las principales 
limitaciones teóricas de este índice de sequía se destaca que no es válido 
cuando el valor de la ETP es igual a 0, lo cual es muy común en regiones frías 
en invierno (Vicente-Serrano et al., 2015). El RDI se expresa en tres formas: 
el valor inicial (a, ), el RDI normalizado (RDIn) y el RDI estandarizado (RDIst). 
El valor inicial («,) se presenta en forma agregada utilizando una escala de 
tiempo mensual, y puede calcularse sobre una base mensual, estacional o 
anual (Asadi-Zarch, Malekinezhad, Mobin, Dastorani, 8 Kousari, 2011), y sus 
valores se calculan mediante la siguiente ecuación: 


0 


(0) 
A 
ds ¡E ETP; 


¡=1,2,33,...,Na (18) 
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donde P;; Y ETP,, son la precipitación y la evapotranspiración potencial del mes 
j del año i, respectivamente, en mm; N,, el número de años de información 
disponible, y k es la escala de análisis en que se calcula el RDI. 


Las formas normalizada y estandarizada del RDI son (Tsakiris, Pangalou, € 
Vangelis, 2007): 


> 
RDIS = 51 (19) 
S, 
(0 yo — Yk 
L k 
RDLico = EE (20) 


donde 3” es la media aritmética de 5” y y” son los logaritmos naturales de 


$; cuya media y desviación estándar son y, y Óyx, respectivamente. 


De acuerdo con diferentes investigaciones con datos de distintos lugares y 
diferentes escalas de tiempo de análisis, se observó que en casi todos los 
casos, los valores de ó, siguen tanto una distribución Gamma como una 
Lognormal. Aunque en la mayoría de los casos la distribución Gamma resultó 
más exitosa (Kousari et a/., 2014; Cai, Zhang, Yao, € Chen, 2015). 


Según Tigkas (2008) y Asadi-Zarch et a/. (2011) se ha probado que el cálculo 
de RDIst podría realizarse mejor ajustando la distribución Gamma a la 
distribución de frecuencias de $,, ya que este método tiende a resolver el 
problema de calcular el RDIst para las escalas de tiempo k pequeñas, que 
pueden incluir valores de precipitación cero (8, = 0), para los cuales la ecuación 
(20) no está definida. 


En este estudio, el RDI se obtuvo ajustando la distribución Gamma a los 
valores 6,, de forma similar al SPI, mediante las ecuaciones (10) a (17), 
sustituyendo los valores de $, en lugar de x. 
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El Índice de Precipitación-Evapotranspiración 
Estandarizado (SPEI) 


La principal limitación del SPI es que se basa sólo en la precipitación e ignora 
otras variables que afectan la demanda de agua atmosférica, como 
temperatura, velocidad del viento, radiación solar y déficit de presión de vapor 
(McEvoy, Huntington, Abatzoglou, € Edwards, 2012). Para superar esta 
limitación, Vicente-Serrano et a/. (2010) desarrollaron el SPEI, el cual combina 
la sensibilidad del PDSI con los cambios en la demanda de la evaporación 
causados por fluctuaciones y tendencias de la temperatura, con la sencillez de 
cálculo y la naturaleza multiescalar del SPI (Vicente-Serrano et al., 2010; 
Banimahd € Khalili, 2013; Li et a/., 2014). 


El SPEI se basa en un balance de agua climático mensual (precipitación menos 
evapotranspiración potencial, llamada serie D) que se calcula en diferentes 
escalas de tiempo y matemáticamente es similar al SPI, pero incluye el papel 
de la evapotranspiración potencial (López-Moreno et al., 2013). Su cálculo 
sigue un enfoque similar al seguido para el cálculo del SPI, pero utilizando una 
distribución de probabilidad Log-logística de tres parámetros en lugar de la 
distribución Gamma de dos parámetros (Figura 2). La clasificación de sequía 
del SPI puede ser usada también para evaluar el SPEI (Tabla 2). 


Para el cálculo del SPEI, primero se calcula la diferencia entre la precipitación 
y la evapotranspiración potencial del mes j del año ¡ (serie DU), en mm; 
después la serie DU se agrega en la escala de tiempo k de interés para obtener 
las series D;?, en milímetros. 


DÚ =P, — ETP;; (21) 


donde P;; y ETP;, tienen el mismo significado que en la ecuación (18). 
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Para el cálculo del SPEI es necesario utilizar una distribución de probabilidad 
de tres parámetros, dado que en una de dos parámetros, la variable x tiene 
un límite inferior de cero (0<x<o0), mientras que en una función de 
distribución de tres parámetros, x puede tomar valores en el rango (y <x<oo, 
donde y es el parámetro de origen de la distribución), por lo cual, Vicente 
Serrano et al. (2010) propusieron utilizar la distribución Log-logística de tres 
parámetros, pues la serie D, puede tener valores negativos, y cuya función de 
densidad de probabilidad es: 


ta] co 


donde a, f, son los parámetros de escala, forma y origen, respectivamente, 
Y 
para los valores D, en el rango de (y > D, < 0). 


Los parámetros de la distribución Log-logística pueden obtenerse por diferente 
métodos, de los cuales, el método de momentos L es el más robusto y de fácil 
aplicación (Ahmad, Sinclair, 8 Werritti, 1988; Singh, Guo, 8i Yu, 1993): 


e ZW; —= Wo 
B ae 6w,; —= Wo — 6w, de 
(wy — 2w1)B 
“TA+1DTA1/0 pe 
1 1 
y =0w ar (147) (1-3) (25) 


donde r(1+1/f) es la función Gamma de (1+ 1/6), y w, es el momento de 
probabilidad pesada de orden s (s = 0,1,2): 
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n 


fa Vs 


i=1 


donde x, es la muestra aleatoria ordenada (x, < x»,--*,< x,) de los valores de 
Dz y n es el tamaño de la muestra. 


La distribución de probabilidad de las series D,, de acuerdo con la distribución 
Log-logística está dada por (Vicente-Serrano et a/., 2010): 


F(x) = E m3 (- ES JN (27) 


Finalmente, para obtener los valores del SPEI, los valores de F(x) se 
transforman en la variable normal estándar Z, sustituyendo F(x) por H(x) en 
las ecuaciones (16) y (17). 


Comparación de los índices de sequía 


Los índices SPI, RDI y SPEI siguen metodologías similares y sus valores tienen 
el mismo significado estadístico, por lo tanto, son comparables. Para realizar 
la comparación entre los índices de sequía se utilizaron los coeficientes de 
correlación de Pearson (r) y de Spearman (r;). 


Resultados y discusión 
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En esta investigación se definió un evento de sequía como aquel en el que los 
valores del SPI, RDI o SPEI toman valores del índice iguales o menores a -1.0. 
Se determinó el número de eventos de sequía, así como la duración y 
magnitud de cada evento para las escalas de 3, 6 y 12 meses, para cada una 
de las nueve estaciones analizadas. 


Un análisis general de los índices de sequía analizados muestra que tanto local 
como regionalmente estos índices siguen el mismo patrón temporal para el 
periodo de 1961-2012. De acuerdo con el SPI, el RDI y el SPEI, los periodos 
húmedos más importantes fueron de 1965-1968, 1990-1992 y 2003-2005, y 
los principales eventos de sequía se produjeron en 1969-1970, 1982-1983, 
1998-2000 y 2010-2012. 


Los tres índices analizados mostraron resultados consistentes en la detección 
de los grandes eventos de sequía respecto a su magnitud y duración (Tabla 
3). Las series del SPI, RDI y SPEI tuvieron pocas diferencias entre sí, 
independientemente de la escala de tiempo de análisis. Estos resultados 
muestran que en condiciones climáticas donde existe baja variabilidad 
interanual de la temperatura, los índices de sequía analizados responden 
principalmente a la variabilidad de la precipitación. 


Tabla 3. Valores máximos de las características de las sequías para todas las 
estaciones y escalas de tiempo. 


Índice | Estación Duración Severidad (magnitud) Intensidad 
Valor Periodo Valor Periodo Valor Periodo 

32052 |6 od “[-8.426 |12/2010 - 05/2011 | -1.483 OS j 

32028 |7 ES “|-9.005 |12/1988 - 06/1989 | -1.361 Pl i 

SPI-3 [32018 6 o “|-9.696 |12/2010 - 05/2011 | -1.713 ió i 

32024 (6 0 “|-9.396 |12/2000 - 05/2001 | -1.574 os i 

32032 (6 Dados “(-8.936 | 12/2010 - 04/2011 | -1.787 RN i 

32052 7 on “|-9.711 | 12/2010 - 06/2011 | -1.483 OS E 
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32028 (7 arta -9,109 |12/1988 - 06/1989 | -1.361 Pe 
RDI-3 [32018 (6 ol -9.877 | 12/2010 - 05/2011 | -1.713 es 
32024 16 tor -9.414 |12/2000 - 05/2001 |-1.574 idol 
32032 16 o -10.102 | 12/2010 - 05/2011 |-1.787 IRON 
32052 |6 a -9.775 [05/1974 - 10/1974 | -1.900 De taoe 
20608. ls dio -8.923 | 05/1980 - 08/1980 | -2.231 Delos 
SPEL-3 (32018 |5 id -7.463 | 05/1962 - 09/1962 | -2.048 de ido 
32024  |4 pts -7.247 [04/1989 - 07/1989 | -1.969 e 
32032 |4 dl -9,134 |04/2011 - 07/2011 |-2.320 la 
32052 10 Dio -14.918 | 12/1984 - 09/1985 |-2.416 de no 
32028 (6 pa -10.971 | 02/1989 - 07/1989 | -2.173 E 
spPI-6 (32018 |5 No -10.564 | 03/2011 - 07/2011 | -2.113 a 
32024 16 edo -12.945 | 01/2001 - 06/2001 |-2.189 DS 00 
32032 |4 a -8.633 |03/2011 - 06/2011 |-2.525 ua 
32052 (15 ne -20.798 | 12/1984 - 02/1986 | -2.437 io 
32028 |5 iodo -10.325 | 03/1989 - 07/1989 | -2.094 De ió 
RDI-6 |32018 (6 eb -12.613 | 03/2011 - 08/2011 | -2.166 pb 
32024 16 e -13.096 | 01/2001 - 06/2001 | -2.243 dai 
32032 |4 oí -9.290 |03/2011 - 06/2011 |-2.622 ae 
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32052 6 losa -10.257 | 06/1974 - 11/1974 | -1.877 Ped 
32028 6 En -12.848 | 06/1980 - 11/1980 | -2.141 on 
SPEI-6 | 32018 7 va -12.414 | 05/2011 - 11/2011 | -1.881 O 
32024 7 oros -12.949 | 04/1989 - 10/1989 | -1.850 O 
32032 9 ee -16.395 | 03/2011 - 11/2011 | -1.822 E 
32052 40 duos -65.884 | 06/1985 - 09/1988 | -1.837 Po 
32028 18 io -30.437 | 07/2011 - 12/2012 | -2.074 sd 
SPI-12 |32018 18 2 -38.158 | 02/2011 - 07/2012 | -2.120 po 
32024 20 tn -35.080 | 08/1989 - 03/1991 | -1.754 en 
32032 23 NOE -44,154 | 02/2011 - 12/2012 | -1.920 re 
32052 40 sio -61.101 | 06/1985 - 09/1988 | -1.832 e 
32028 14 dao -24.325 | 07/2011 - 08/2012 | -1.835 o 
RDI-12 | 32018 18 aa -36.955 | 02/2011 - 07/2012 | -2.053 ne 
32024 20 Ea -35.433 | 08/1989 - 03/1991 | -1.772 bea 
32032 23 E -52.654 | 02/2011 - 12/2012 | -2.289 no 
32052 25 oe -40.618 | 12/1961 - 12/1963 | -1.625 o 
32028 12 peo -22.007 | 06/1980 - 05/1981 | -1.935 es 
SPEI-12 | 32018 18 o -30.817 | 02/2011 - 07/2012 | -1.814 o 
32024 24 o -38.493 | 06/1989 - 05/1991 | -1.629 e 
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11/2010 - 11/2010 - 
32032 26 12/2012 50.699 | 11/2010 - 12/2012 | -1.950 12/2012 


Los resultados de los índices SPI y RDI son similares en las escalas de 3 y 6 
meses (SPI-3 y RDI-3; y SPI-6 y RDI-6). Con el SPI-3 y el RDI-3 la mayoría 
de los periodos de sequía se detectó en los meses de enero a mayo; con el 
SPEI-3 se detectó de abril a julio. Respecto al SPI-6 y RDI-6, la mayoría de los 
periodos secos ocurrió en los meses de marzo a junio; con el SPEI-6 se 
presentó de mayo a septiembre. Se observa que, para estas escalas, el SPEI 
detecta la sequía varios meses después que el SPI y el RDI debido a la relación 
que existe entre precipitación y temperatura en los primeros meses del año. 
En la escala de 12 meses, los resultados de los tres índices son similares, 
aunque en los periodos importantes en cuanto a duración, magnitud e 
intensidad presentan mayor semejanza el SPI-12 y el RDI-12. Según el SPI- 
12, el RDI-12 y el SPEI-12 en 1962, 1969, 1979, 1982, 1998, 1999, 2000, 
2011 y 2012 hubo importantes periodos de sequía. Años que en su mayoría 
fueron detectados con sequía por Campos-Aranda (2014) en la estación 
climatológica Fresnillo en el centro de Zacatecas, con el Índice del Déficit 
Probabilístico de Precipitación, el SPI y el RDI. Campos-Aranda (2016) también 
reportó sequías en una escala anual en gran parte de Zacatecas con el índice 
de aridez en estos mismos años. El periodo de sequía más importante en 
cuanto a intensidad se presentó entre 2011 y 2012. 


La evolución de las series SPI, RDI y SPEI en la escala de tres meses (corto 
plazo) presentó una frecuencia temporal alta de los periodos de sequía y de 
humedad. Conforme la escala de tiempo se incrementa a 6 y 12 meses, los 
periodos secos y húmedos mostraron una frecuencia temporal menor y una 
mayor duración (Figura 3). Se definen con mayor claridad las temporadas 
secas y húmedas, sobre todo en la escala de 12 meses. 
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Figura 3. Evolución temporal de la sequía en las escalas de 3, 6 y 12 meses 
en la estación El Sauz (32018). 


Características medias de las sequías 


En la Tabla 4 se presenta un resumen del número de eventos de sequía y su 
duración media, para periodos de sequía a corto, mediano y largo plazos, 
determinados con los tres índices de sequía evaluados. El número medio de 
episodios de sequía disminuyó con el aumento de las escalas de tiempo. 


Tabla 4. Número de eventos y (duración media) de los eventos de sequía 
determinados por el SPI, el RDI y el SPEI para el periodo 1961-2012. * El 
registro de esta estación es de 1971-2012. 


SPEI- SPEI- SPI- RDI- SPEI- 


Estación | SPI-3 | RDI-3 3 SPI-6 | RDI-6 6 12 12 12 


74, 


pa 43 44 38 28 30 27 16 14 19 
Q4 1Q2 162 (60D (60 ((41) |(68) |(81) |(6.1) 
aa 36 41 47 41 39 34 21 16 19 
(28) 1(QN 1(Q3) 1Q60 1Q) 1(Q9 |(51) |(58) |(5.1) 
20 37 42 45 32 26 34 10 11 19 
(2.9) 1(2.5) 1(2.6) 1(3.0) 1(3.2) ((3.2) |(6.2) |(5.6) |(4.6) 
204 42 46 47 38 34 36 13 12 13 
(2.3) 1(Q3 (QDD (GD (G0 (161) |(89) |(9.4) |(9.1) 
2 48 49 51 37 36 35 18 19 18 
(2.3) 1Q4 1(Q2 1Q9 102.8) 1(3.0) |(5.7) |(51) |(5.8) 
0 41 41 48 39 33 33 19 20 21 
(2.3) 1(2.3) ((Q.0) 1(2.8) 1(3.0) 1(3.3) |(5.2) |(5.0) |(5.4) 
a al 34 40 26 26 28 16 19 18 
(28) 1(Q5) (Q0D 16D 1(G0 1633) 1654) |(48) |(5.3) 
a 47 46 44 38 37 39 18 15 13 
(2.5) (6) (Q1D) 1Q8 1Q9 12.9) |(5.7) |(6.7) |(8.2) 
a 44 47 52 34 31 37 19 21 18 
(2.5) 1(Q49 1(Q0 160 (6D (GD 148 142 |(5.7) 


El número de eventos de sequía promedio en una escala de tiempo a corto 
plazo varió entre 36 y 48 para el SPI-3, de 41 a 49 para el RDI-3 y de 38 a 52 
para el SPEI-3. La duración media de los eventos de sequía osciló de 2.3 a 
2.9, de 2.2 a 2.7 y de 2.0 a 3.2 meses para el SPI-3, el RDI-3 y el SPEI-3, 
respectivamente. En 70% de las estaciones climatológicas el evento de sequía 
de mayor magnitud alcanzó la categoría de sequía severa, con el SPI-3 y el 
RDI-3, mientras que con el SPEI-3, 90% de las estaciones climatológicas 
alcanzó la categoría de sequía extrema. La duración máxima de los periodos 
de sequía registrados con el SPI-3 y el RDI-3 varió de 5 a 8 meses, mientras 
que para el SPEI-3 dicha duración fue de 4 a 7 meses. 


En las dos últimas décadas se presentaron los eventos de sequía más intensos 
en la escala de tres meses. Según el SPI-3 y el RDI-3, estos eventos fueron 
de febrero a mayo de 1999, y de enero a mayo de 2011, y de acuerdo con el 
SPEI-3 éstos fueron de mayo a julio de 2011 y de mayo a junio de 1998. 
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Para la escala de seis meses, los eventos de sequía detectados por los tres 
índices presentan resultados relativamente similares. Los eventos de sequía 
promedio variaron de 28 a 41, de 26 a 39 y de 27 a 39 para el SPI-6, el RDI- 
6 y el SPEI-6, respectivamente. La duración media de tales eventos va desde 
2.6 a 3.1 meses para el SPI-6, de 2.7 a 3.2 meses para el RDI-6 y de 2.9 a 
4,1 meses para el SPEI-6. Los eventos de mayor intensidad alcanzaron la 
categoría de sequía extrema con el SPI-6 y el RDI-6 en todas las estaciones; 
con el SPEI-6 sólo 20% alcanzó esta categoría, las demás estaciones 
presentaron sequía severa. De acuerdo con el SPEI-6, las duraciones máximas 
de los periodos secos fueron muy similares entre estaciones, variando de 6 a 
9 meses, mientras que de acuerdo con el SPI-6 y el RDI-6 presentaron mayor 
variación, de 4 a 15 meses. 


Con base en el SPI-6, el RDI-6 y el SPEI-6, la región noreste de Zacatecas 
presentó sequías importantes en cuanto a magnitud e intensidad en las 
décadas de 1970 y 1980. Para las décadas de 1990 y 2000, las sequías 
importantes se presentaron en el centro y sur del estado. 


En la escala de seis meses, el RDI-6 detectó el periodo de sequía de mayor 
duración (15 meses), el cual ocurrió en la región norte de Zacatecas (estación 
32052). El periodo de sequía más largo detectado por el SPI-6 fue en la 
estación 32001 (norte de Zacatecas), también detectado por el RDI-6, el cual 
duró 13 meses, de junio de 1987 a junio de 1988. El SPEI-6 detectó sequía 
para este mismo periodo y estación, pero en dos subperiodos, uno de seis 
meses (de julio a diciembre de 1987) y otro de tres (de mayo a julio de 1988). 
El periodo de sequía de mayor duración detectado con el SPEI-6 fue de nueve 
meses en las estaciones 32032 y 32054, de marzo a noviembre de 2011, sur 
y centro-oeste de Zacatecas, respectivamente. En este mismo periodo, el SPI- 
6 y el RDI-6 detectaron un periodo de sequía de sólo cuatro meses, de marzo 
a junio de 2011. 


Respecto a la escala de tiempo de 12 meses (mediano y largo plazo), el 
número de eventos de sequía promedio es similar para los tres índices, 
variando de 10 a 21 para el SPI-12, de 11 a 21 para el RDI-12 y de 13 a 21 
para el SPEI-12. La duración media de los eventos de sequía oscila de 4.8 a 
8.9 para el SPI-12, de 4.2 a 9.4 para el RDI-12 y de 4.6 a 9.1 meses para el 
SPEI-12. En esta escala, los eventos de mayor intensidad alcanzaron la 
categoría de sequía extrema en 70% de las estaciones con el SPI-12 y el RDI- 
12, y 30% de sequía severa, mientras que con el SPEI-12, sólo el 10% de las 
estaciones presentó sequía extrema, el 90% restante presentó sequía severa. 
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Las duraciones máximas de los periodos secos mediante el SPI-12 y el RDI-12 
fueron similares, variando de 16 a 40 y de 14 a 40 meses, respectivamente; 
estas duraciones fueron menores para el SPEI-12, variando de 12 a 27 meses. 


Durante el periodo de 52 años de registros climatológicos analizados en el 
estado de Zacatecas, las intensidades máximas de sequía en la escala de 
tiempo anual fueron detectadas por el SPI, lo cual coincide con los resultados 
obtenidos por Campos-Aranda (2014, 2015) al analizar las sequías con el SPI 
y el RDI en el centro de Zacatecas; en las escalas de 3 y 6 meses, las 
intensidades máximas de sequía fueron detectadas por el SPEI y el RDI, 
respectivamente. 


Las sequías de mayor intensidad para las tres escalas de tiempo se 
presentaron entre 1998 y 2011. El RDI y el SPEI también identificaron sequías 
extremas en 1974, 1979, 1982 y 1989. En la Figura 4 se muestra la 
distribución espacial del SPI, el RDI y el SPEI en el mes de mayor sequía en 
las tres escalas de tiempo analizadas. En la escala de tres meses, el SPI-3 y 
el RDI-3 identificaron el mes de abril de 2011 como el de mayor intensidad de 
sequía; dichos índices tienen un comportamiento espacial prácticamente igual 
(Figura 4a y Figura 4b). La mayor parte de Zacatecas (82%) presentó sequía 
severa y sólo el 8% sequía extrema. En junio de 2011 se presentó la sequía 
de mayor intensidad con el SPEI-3, con sequía extrema en 59% del estado, 
principalmente en la parte centro (Figura 4c). En la escala de seis meses, abril 
de 2011 también fue el mes de mayor intensidad de sequía con el SPI y RDI, 
el cual presentó sequía extrema en prácticamente todo el territorio de 
Zacatecas (figura 4d y Figura 4e). Julio de 2011 fue el mes de mayor sequía 
identificado por el SPEI-6; la mayor parte del territorio presentó sequía 
extrema (49%) y severa (39%). La sequía extrema ocurrió sobre todo en el 
centro y noreste de Zacatecas. En la escala de 12 meses, septiembre de 2011 
fue el mes de mayor intensidad de sequía identificado por los tres índices. El 
SPI-12 y el RDI-12 tuvieron un comportamiento espacial semejante: 
identifican alrededor de 55% de Zacatecas con sequía extrema (Figura 4g y 
Figura 4h), principalmente sur y oeste del estado. El SPEI-12 identificó sequía 
extrema en 33% de Zacatecas, en las mismas regiones que el SPI-12 y el RDI- 
12, pero con una extensión menor. 
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Figura 4. Distribución espacial de los índices SPI, RDI y SPEI en el mes de 
sequía más intensa para las escalas de 3, 6 y 12 meses. 


Frecuencia de las sequías 


La frecuencia de las sequías se analizó con base en las diferentes clases de 
sequía (Tabla 2) relacionadas con los valores de los índices de sequía (SPI, 
RDI y SPEI). En la escala de tres meses, para los índices de sequía analizados, 
el tiempo promedio (porcentaje de frecuencia promedio) de periodos secos y 
húmedos fue de 17.1 y 16.9%, respectivamente. El 66.0% del tiempo presentó 
condiciones normales de precipitación. El SPI-3 y el RDI-3 presentaron mayor 
frecuencia de sequías moderadas (12.5 y 12.2%, respectivamente) que el 
SPEI-3, mientras que la frecuencia de sequías severas fue menor para el SPI- 
3 y el RDI-3 (4.1 y 4.8%, respectivamente) que para el SPEI-3. De los tres 
índices, el SPI-3 y el RDI-3 presentaron el menor número de sequías extremas. 
En la Figura 5a se observa el porcentaje de frecuencia de periodos secos y 
húmedos para el SPI-3, el RDI-3 y el SPEI-3. Con el SPI-3 y el RDI-3 sólo la 
región centro-oeste de estudio (estación 32054) presentó sequías extremas 
(14 y 15 eventos, respectivamente), mientras que con el SPEI-3 todas las 
estaciones presentaron eventos de sequía extrema, variando de 3 a 16 
eventos para las estaciones 32052 y 32032, norte y sur de Zacatecas, 
respectivamente. 
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Figura 5. Distribución de frecuencia de los valores de los índices SPI, RDI y 
SPEI: a) escala de tres meses; b) escala de seis meses; y c) escala de 12 
meses, para el periodo de 1961-2012. 


En la escala de seis meses, según los tres índices multiescalares, en promedio 
67.4% del tiempo hubo condiciones normales de precipitación. La frecuencia 
de los eventos secos fue de 16.7%, y de los eventos húmedos de 15.9%. El 
SPEI-6 presentó mayor frecuencia de sequías moderadas (11.0%) que el SPI- 
6 y el RDI-6 (8.7 y 7.7%, respectivamente). Lo mismo ocurrió con las sequías 
severas (Figura 5b). Sin embargo, el SPI-6 y el RDI-6 presentaron mayor 
frecuencia de sequías extremas que el SPEI-6, 3.4 y 3.8%, respectivamente. 
Con el SPI-6 y el RDI-6 se detectaron 26 y 29 eventos de sequía extrema, los 
cuales ocurrieron en el centro de Zacatecas (estaciones 32018 y 32030, 
respectivamente), mientras que con el SPEI-6, el número máximo de eventos 
de sequía extrema fue de 15, aconteciendo en la región norte (estación 
32001). Por último, en la escala de 12 meses, el porcentaje de frecuencia 
promedio, de eventos normales, húmedos y secos fue de 66.6, 16.9 y 16.5%, 
respectivamente. La frecuencia de eventos de sequía para el SPI-12, el RDI- 
12 y el SPEI-12 fueron 16.3, 15.8 y 17.4%, respectivamente. El 
comportamiento y los valores de la frecuencia de las diferentes clases de 
sequía para el SPI-12, el RDI-12 y el SPEI-12 son similares a los que presentan 
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el SPI-6, el RDI-6 y el SPEI-6. Es decir, el SPEI-12 presenta mayor frecuencia 
de ocurrencia de sequías moderadas y severas que el SPI-12 y el RDI-12, y 
para sequías extremas sucede lo contrario, como se observa en la Figura 5c. 
Con el SPI-12 y el RDI-12 se detectaron 23 y 22 eventos de sequía extrema 
en la estación 32032 (sur de Zacatecas), mientras que con el SPEI-12, el 
número máximo de eventos de sequía extrema fue de 17 en la estación 32001 
(norte de Zacatecas). 


Correlación entre los índices de sequía 


Considerando las tres escalas de tiempo y las nueve estaciones analizadas, los 
índices SPI vs. RDI son los índices que tienen correlación más alta, y sus 
coeficientes de correlación r varían de 0.920 a 0.996. Los coeficientes de 
correlación más bajos los presentan los índices SPI vs. SPEI, variando de 0.336 
a 0.982. En la Figura 6 se muestra la tendencia general de las relaciones entre 
los índices SPI vs. RDI, SPI vs. SPEI y RDI vs. SPEI, en las escalas de 3, 6 y 
12 meses mediante diagramas de dispersión para las series de los índices de 
la estación El Sauz (32018), ubicada en la región centro de Zacatecas. 


81 


Figura 6. Diagramas de dispersión para las series de SPI vs. RDI, SPI vs. 
SPEI y RDI vs. SPEI para la estación El Sauz (32018), para el periodo de 
1961-2012. 


En la escala de tres meses, los coeficientes de correlación r entre SPI-3 vs. 
RDI-3 varían de 0.920 a 0.969, mientras que entre SPI-3 vs. SPEI-3 y RDI-3 
vs. SPEI-3 varían de 0.336 a 0.804 y de 0.583 a 0.867, respectivamente. En 
general, las correlaciones para la escala de seis meses fueron ligeramente más 
altas que para la escala de tres meses, variando de 0.943 a 0.964 para SPI-6 
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vs. RDI-6, de 0.410 a 0.834 para SPI-6 vs. SPEI-6, y de 0.643 a 0.906 para 
RDI-6 vs. SPEI-6. En la escala de 12 meses, los coeficientes de correlación 
varían de 0.952 a 0.996 para SPI-12 vs. RDI-12, de 0.926 a 0.982 para SPI- 
12 vs. SPEI-12, y de 0.952 a 0.993 para RDI-12 vs. SPEI-12. Se observó que 
para las escalas y estaciones analizadas, a menor escala, menor correlación 
entre los índices (Tabla 5). En las escalas de 3 y 6 meses, las correlaciones 
más bajas se presentaron al norte del estado de Zacatecas (estación 32052) 
y las más altas al sur (estación 32032). 


Tabla 5. Coeficientes de correlación r y r, entre los índices de sequía SPI, 


RDI y SPEI. 
Esta- 
ción 32001 32003 32028 32018 32032 32054 32052 32030 32024 
Coefi- 
Índice | ciente r T r T r T ñ T r T r Ts r Ts r Te r TE 
SPI-3 RDI-3 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 
33 28 55 48 48 50 47 46 69 69 49 48 20 26 34 40 59 56 
SPEI-3 0.4 0.3 0.6 0.6 0.5 0.5 0.6 0.6 0.8 0.8 0.7 0.7 0.3 0.2 0.7 0.7 0.6 0.6 
65 82 71 75 62 23 36 33 04 26 65 ZO 36 41 66 73 33 22 
RDI-3 | SPEI-3 0.6 0.6 0.7 0.8 0.7 0.7 0.7 0.8 0.8 0.9 0.8 0.9 0.5 0.5 0.8 0.8 0.7 0.7 
84 54 93 34 25 17 78 11 67 10 62 06 83 31 62 99 66 83 
SPI-6 RDI-6 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 
43 33 58 31 47 35 49 29 62 39 52 36 44 41 43 11 64 41 
SPEI-6 0.4 0.4 0.6 0.7 0.5 0.5 0.6 0.6 0.8 0.8 0.8 0.8 0.4 0.3 0.7 0.8 0.6 0.6 
99 45 83 05 60 68 37 55 34 75 o1 21 10 31 71 00 54 77 
RDI-6 | SPEI-6 0.7 0.7 0.8 0.8 0.7 0.7 0.8 0.8 0.9 0.9 0.9 0.9 0.6 0.6 0.8 0.9 0.7 0.8 
16 04 12 91 50 82 03 58 06 70 oi 45 43 00 86 48 90 60 
SPI- RDI-12 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 
12 96 96 96 96 87 88 96 95 52 36 91 90 96 96 80 77 91 89 
SPEI-12 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 0.9 
64 80 82 86 50 51 81 85 26 12 76 81 27 40 52 59 62 60 
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0.9 [09 [0.9 [0.99 [0.9 [0.9 [0.9 |0.9 [0.9 |0.9 [0.9 [0.9 [0.9 [0.9 [0.9 | 0.9 
78 92 93 97 86 85 93 97 92 97 91 98 52 60 91 97 
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0.9 
89 


Los valores de los coeficientes de correlación de Spearman (Tabla 5) 
presentaron valores semejantes a los del coeficiente de Pearson. Tal como 
sucedió con el coeficiente r, se observó que para las escalas analizadas los 
coeficientes r, aumentaban con la escala. También se observó que para las 
escalas de 3 y 6 meses, los índices SPI vs. RDI presentaron los mayores 
valores de r;, en tanto que los índices SPI vs. SPEI presentaron los menores 
valores. Para la escala de 12 meses, los índices RDI-12 vs. SPEI-12 y SPI-12 
vs. RDI-12 presentaron los mayores coeficientes r,. Los coeficientes r y 1; 
presentaron mayor similitud conforme aumentaba la escala. 


Para el caso de estudio, el empleo del SPI y del RDI resultó más adecuado 
para la caracterización de la sequía meteorológica. Según Gocic y Trajkovic 
(2014), y Xu et al. (2015), el RDI y el SPI son índices más adecuados que el 
SPEI para caracterizar la sequía en zonas áridas y semiáridas, lo cual coincide 
con los resultados de esta investigación. 


Por otro lado, cabe destacar que la selección de la función de distribución de 
probabilidad en el análisis de los índices de sequías, puede influir en los 
resultados obtenidos, subestimando o sobreestimando los eventos de sequía. 
Por lo tanto, es conveniente determinar mediante un análisis de frecuencia, la 
función de distribución que se ajuste mejor a las series de precipitación, 
precipitación/evapotranspiración y precipitación - evapotranspiración, para 
evaluar los índices SPI, RDI y SPEI, respectivamente. 


Conclusiones 


Los índices multiescalares de sequía meteorológica utilizados en esta 
investigación permitieron caracterizar las sequías en los últimos 50 años en el 
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estado de Zacatecas, México, obteniendo el número, inicio, final, duración, 
magnitud, intensidad y categoría de los eventos de sequía. 


Los tres índices de sequía aplicados mostraron resultados consistentes en la 
detección de los eventos de sequía, en particular los de mayor magnitud; con 
independencia de la escala temporal del índice, su uso permitió un análisis 
complementario de la variabilidad espacial y temporal de las sequías. 


El número de eventos de sequía detectados con los tres índices tuvo menor 
frecuencia y mayor duración conforme la escala temporal de análisis aumentó. 
Las tres escalas presentaron eventos de sequía de categoría moderada a 
extrema, de corta y mediana duración en las escalas de 3 y 6 meses, y de 
media a larga duración para la escala de 12 meses. El RDI presentó el mayor 
número de sequías en la escala de tres meses, y en la escala de 6 y 12 meses, 
el SPI y el SPEI, respectivamente. 


El análisis espacial de las sequías muestra que las áreas más afectadas por 
este fenómeno son la región noroeste, centro y sur del estado. En estas zonas 
se encuentran ubicadas las principales zonas agrícolas de temporal y de riego 
de Zacatecas. 


Los coeficientes de correlación de Pearson y de Spearman fueron mayores 
conforme aumentó la escala de análisis de los índices multiescalares. Los 
coeficientes r y r, confirman que existe una mayor similitud entre los 
resultados del SPI y el RDI para las escalas de 3 y 6 meses. En la escala de 
12 meses, el comportamiento de los tres índices fue semejante, siendo los 
resultados del RDI y el SPEI los que presentaron la mayor similitud. Por lo 
cual, si se dispone de información de temperatura es conveniente aplicar el 
RDI para el análisis de sequías en regiones semiáridas, pues incorpora otras 
variables climatológicas que juegan un papel importante durante el desarrollo 
de un evento de sequía, como la evapotranspiración. No obstante, en ausencia 
de información de temperatura o por facilidad de cálculo, para las escalas de 
análisis el SPI es una buena alternativa para la estimación de las 
características de las sequías. 


La caracterización de las sequías históricas es importante para generar 
información relacionada con déficits hídricos durante los periodos de sequía, 
la cual permite diseñar políticas con el fin de disminuir los daños de eventos 
futuros de sequía. Dada la importancia del sector agrícola en el estado de 
Zacatecas, los resultados encontrados pueden ayudar principalmente en la 
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toma de decisiones sobre los patrones de cultivo y la asignación óptima de 
agua en el Distrito de Riego 034 Estado de Zacatecas. 
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